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COEULERIAN GRAPHS 


MATTHEW FARRELL AND LIONEL LEVINE 


Abstract. We suggest a measure of “Eulerianness” of a finite directed graph and 
define a class of “coEulerian” graphs. These are the graphs whose Laplacian lattice 
is as large as possible. As an application, we address a question in chip-firing posed 
by Bjorner, Lovasz, and Shor in 1991, who asked for “a characterization of those 
digraphs and initial chip configurations that guarantee finite termination Bjorner 
and Lovasz gave an exponential time algorithm in 1992. We show that this can 
be improved to linear time if the graph is coEulerian, and that the problem is 
NP-complete for general directed multigraphs. 


1. Introduction 

In this paper G = (V, E) will always denote a finite directed graph, with loops and 
multiple edges permitted. We assume throughout that G is strongly connected: for 
each v, w G V there are directed paths from v to w and from w to v. Trung Van 
Pham [22] introduced the quantity 

Mq = gcd{/t(?;)|u G V} 

where n(v) is the number of spanning trees of G oriented toward v. We will see that 
Mq, which we will call the Pham index of the graph G, can be interpreted as a 
measure of “Eulerianness”. 

A finite directed multigraph G is called Eulerian if it has an Eulerlian tour (a 
closed path that traverses each directed edge exactly once). We are going to take the 
view that Eulerianness is an algebraic property of the graph Laplacian A acting 
on integer-valued functions / G Z v by 

A/0) = d v f[y) - X /(tail(e)). (1) 

head(e)=t; 

Here d v is the outdegree of vertex v. The context is the following well-known equiv¬ 
alence, where 1 denotes the constant function l(u) = 1 for all v G V. 
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Proposition 1.1. The following are equivalent for a strongly connected directed 
multigraph G = (V,E). 

(1) ker(A : Z v Z v ) = Zl. 

(2) Mq = k(s) for all s € V. 

(3) G is Eulerian. 

Our main result is in some sense dual to Proposition 1.1: it gives several equivalent 
characterizations of the graphs with Pham index 1. These coEulerian graphs are 
the farthest from being Eulerian. 

Our motivation for considering coEulerian graphs and the Pham index comes from 
chip-firing, which we now describe. A chip configuration on G, or simply configu¬ 
ration for short, is a function a : V —>• Z. If cr(v) > 0 we think of a pile of <r(v) chips 
at vertex v, and if cr(v) < 0 we think of a hole waiting to be filled by chips. Denoting 
by d v the outdegree of vertex v, we say that v is stable for a if cr(v) < d v , and active 
for a otherwise. A vertex v can fire by sending one chip along each outgoing edge, 
resulting in the new configuration 

o' = a — A b, 0 

where A is the graph Laplacian (1) and 8 v (w ) is 1 if v = w and 0 otherwise. Con¬ 
cretely, we may think of a, a' as column vectors and A 5 V as a column of the matrix 

» _ j d wv , v w 

^VW — \ , , 

[dv-dwy V = W 

where d wv denotes the number of directed edges of G from w to v. More generally, 
we can specify a firing vector x € N' and fire each vertex v a total of x(u) times, 
resulting in the configuration o' = a — Ax. Here and throughout, N = {0,1,2,...}. 

A legal firing sequence is a finite sequence of configurations <7o, - - -, cr^ such that 
each ori for i = 1 ,k is obtained from cq_i by firing a vertex that is active for 
ai —i . A configuration a is called stable if a(v) < d v for all v E V. We say that 
a stabilizes if there is legal firing sequence a = <ro,..., 07 - such that cr/ ; is stable. 
Bjorner, Lovasz, and Shor posed the following problem in 1991 [ ]. 

The halting problem for chip-firing: 

Given the adjacency matrix of a finite, strongly connected multigraph G 
and a chip configuration a on G with a > 0, 

Decide whether a stabilizes. 

Write \<t\ = cr(v) for the total number of chips. This quantity is conserved 

by firing (since |A<5„| = 0 for all v € V). The maximal stable configuration 

^max(^) = d v 1 

has |Umax| = — ffV. By the pigeonhole principle, any configuration a with 

|oj > jfE — #V has at least one unstable vertex, so such a does not stabilize. A 
natural question arises: Which directed graphs have the property that every chip 
configuration of JfE — ffV chips stabilizes? These graphs are the subject of our main 
result. 
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We write Zq for the set of cr € Z) such that |cr| =0. 

Theorem 1.2. The following are equivalent for a strongly connected directed multi¬ 
graph G = (V, E). 

(1) Im(A : -+ Z v ) = 2%. 

(2) M g = 1. 

(3) A chip configuration a on G stabilizes if and only if |cr| < #E — ffV. 

(4) For all s £ V, the sandpile group K(G,s ) is cyclic with generator f3 s . 

(5) For some s £ V, the sandpile group K(G,s ) is cyclic with generator f5 s . 

Items (1) and (2) are in some sense dual to their counterparts in Proposition 1.1, 
so we propose the term coEulerian for a graph satisfying the equivalent conditions 
of Theorem 1.2. The sandpile group K(G,s) and /3 S are defined below in Section 2. 
For a dual counterpart to item (4), see Proposition 2.12(4). 
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Figure 1. Example of a coEulerian graph: a path of length n with 
edge multiplicities 2 to the right and 3 to the left. It has n(y) = 
tt(v) = 2 v 3 n ~ v spanning trees oriented toward v, so its Pham index is 
M g = gcd(2 n , 2 n_1 3,..., 3 n ) = 1. 


1.1. History. Let us call a graph bidirected if it is obtained from an undirected graph 
by replacing each undirected edge {u,v} by a pair of directed edges (u,v) and (y,u). 
All bidirected graphs are Eulerian. 

One of the earliest results in chip-firing is the following observation of Tardos. 

Lemma 1.3. [26, Lemma 4] Let a be a configuration on a bidirected graph G. If 
there is a legal firing sequence for a in which every vertex of G fires at least once, 
then o does not stabilize. 

Tardos used Lemma 1.3 to prove that for any configuration a on a simple bidirected 
n-vertex graph, if a stabilizes then it does so in 0(n 4 ) firings. Eriksson showed, 
however, that on a general directed graph a configuration may require an exponential 
number of firings to stabilize [13]. Bjorner and Lovasz [5] generalized the “at least 
once” condition of Lemma 1.3 to directed graphs as follows. 

Lemma 1.4. [ ] For every strongly connected multigraph G there is a unique primi¬ 
tive it £ such that A 7 r = 0 .If there is a legal firing sequence for a in which every 
vertex v fires at least ir(v) times, then a does not stabilize. 

This gives a procedure for deciding the halting problem for chip-firing: per¬ 
form legal firings in any order until either you reach a stable configuration or the 
criterion of Lemma 1.4 certifies that a will not stabilize. There is only one problem: 
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the values tv(v) may be exponentially large. Figure 1 shows a coEulerian graph on 
vertex set {0,1,..., n} with tt(v) = 2 v 3 n ~ v . The algorithm just described would run 
for exponential time on this graph, but Theorem 1.2 gives a much faster algorithm to 
decide the halting problem for chip-firing on any coEulerian graph: count the total 
number of chips and compare to ffE — ffV. As far as we are aware, this is the first 
progress on the halting problem for chip-firing on directed graphs since the work of 
Bjorner and Lovasz [ ]. 

1.2. Related work. Pham [ ] introduced the index Mq to answer a question posed 

in [17]: Which directed graphs G have the property that all unicycles of G lie in the 
same orbit of the rotor-router operation? He showed that G has this property if and 
only if Mq = 1, and that in general the number of orbits is Mq. 

The halting problem FOR chip-firing is a special case of the halting problem 
for a class of automata networks called abelian networks. A polynomial time algorithm 
to decide if a given abelian network halts on all inputs appears in [8], where it is 
remarked that the problem of deciding whether a given abelian network halts on a 
given input is a subtler problem. The halting problem FOR chip-firing is of this 
latter type (the “input” to the abelian network is the chip configuration a). 

1.3. Outline. The next section is devoted to the proofs of Proposition 1.1 and Theo¬ 
rem 1.2, and concludes with Proposition 2.13 characterizing the graphs that are both 
Eulerian and coEulerian. In Section 3 we show that despite its being easy for Euler- 
ian graphs and coEulerian multigraphs, the halting problem for chip-firing on 
finite directed multigraphs is NP-complete in general. One ingredient in the proof is 
Theorem 3.1, which expresses an arbitrary (n — l)-dimensional lattice in Zq as the 
Laplacian lattice of a strongly connected multigraph. 

2. Sandpiles and the halting problem 

To prove Theorem 1.2 we will compare chip-firing with and without a sink vertex. 
This kind of comparison appears also in the study of the abelian sandpile threshold 
state [19], and in the extension of the Biggs-Merino polynomial to Eulerian graphs 
[ ] and to all strongly connected graphs [10]. Sections 2.1 and 2.2 review the rele¬ 

vant background on chip-firing and the sandpile group. In Section 2.3 we relate the 
sandpile groups with and without sink, and in Section 2.4 we prove the results stated 
in the introduction. 

2.1. Background. The following result frees us from considering only legal firing 
sequences in looking for an answer to the halting problem for chip-firing. 

Lemma 2.1. (Least Action Principle, [ , Lemma 4.3]) Let a be a chip configuration 
on a finite directed graph. Then a stabilizes if and only if there exists an x £ N v 
such that a — Ax is stable. 

A sizable portion of the ground soon to be covered is motivated by the following 
principle: in looking for a stabilizing firing sequence, instead of firing willy-nilly we 
can establish some structure by designating a special vertex s as the sink, which fires 
only if no other vertex is active. We fire active, nonsink vertices until all nonsink 
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vertices are stable. At this point if the sink is stable we are done; otherwise, we fire 
the sink (once) and repeat. 

The reduced Laplacian A s is the matrix obtained by deleting the row and column 
of A corresponding to the sink s. To emphasize the distinction between A and A s , we 
will sometimes refer to A as the total Laplacian. In what follows we will sometimes 
identify the vertex set V with {1,... , n} and set s = n. 

Definition 2.2. Let G = iV : E) be a finite strongly connected multigraph and fix 
s£k. The sandpile group of G with sink s is the group quotient 

K(G,s) = Z n - 1 /A s Z n - 1 
where A S TA~ X is the integer column-span of A s . 

A sandpile is a chip configuration i] £ Z^ 1 indexed by the nonsink vertices. 
When we wish to emphasize that a chip configuration is defined also at the sink, we 
call it a total configuration. One can imagine that a sandpile is composed of sand 
grains which behave just like chips except that they are small enough to disappear 
down the sink. The definitions “stable” and “firing vector” have obvious analogues 
for sandpiles: a sandpile ij is stable if r/(i) < di for all u* ^ s; and firing vectors 
for sandpiles live in Z n . The sandpile group treats two sandpiles as equivalent if 
one can be obtained from the other by firing nonsink vertices. We write fj for the 
equivalence class of r/ in I\(G , s ). 

On a strongly connected graph, every sandpile stabilizes, and its stabilization does 
not depend on the order of firings [17, Lemmas 2.2 and 2.4]; we denote the stabi¬ 
lization of r] by rf. Next we recall the connection between sandpiles and spanning 
trees. 

Definition 2.3. An oriented spanning tree of a directed graph G = (V, E ) rooted 
at s £ V is a spanning subgraph T = (V, A) such that 

(1) Every vertex v/s has outdegree 1 in T. 

(2) s has outdegree 0 in T. 

(3) T has no oriented cycles. 

Hence an oriented spanning tree has as its limbs edges that point toward the root. 
Let k(s) denote the number of oriented spanning trees in G rooted at s. 

Theorem 2.4. (Matrix tree theorem [25, Theorem 5.6.8] and [ , Lemma 2.8]) For 

a finite strongly connected multigraph G and a vertex s, 

k(s) = det A s = #K(G,s). 

Note that if G is strongly connected then it has at least one spanning tree rooted 
at s, so A s is invertible; since the rows of A sum to 0, this implies that A has rank 
n — 1. 

There is a natural representative for each equivalence class of K(G, s ). To describe 
this representative, we say that a sandpile rj is accessible if from any other sandpile 
it is possible to obtain r/ by adding a nonnegative number of sand grains at each 
vertex and then selectively firing active vertices. A sandpile that is both stable and 
accessible is called recurrent. 
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Theorem 2.5. [17, Cor. 2.16] The set Rec(G, s) of all recurrent sandpiles is an 
abelian group under the operation 

:= {V + C)° 

and it is isomorphic via the inclusion map to the sandpile group I\(G,s). 

The recurrent identity element e s £ Rec(G, s ) is the unique recurrent sandpile 
in A s Z n_1 . The recurrent representative r] vec of a sandpile r) can be found by adding 
the identity and stabilizing: 

Vrec = {v + e s )°. 

Dhar’s burning test [ I ] determines whether a given sandpile on an Eulerian graph 
is recurrent. Speer [ ] generalized the burning test to directed graphs. Dhar’s and 

Speer’s tests are closely related to Lemmas 1.3 and 1.4 respectively. 

2.2. Cyclic subgroups of the sandpile group. For s,v £ V let /3 s (v) = d sv , the 
number of directed edges from s to r. In accordance with our principle of controlled 
sink firing, given a recurrent sandpile ry we are interested in 

C v = {(r, + k(d s )° : k £ N}, 

the set of sandpiles obtainable from 77 by firing the sink s some nonnegative number 
of times and then stabilizing. Note that starting with a recurrent sandpile, adding 
sand grains to the nonsink vertices and then stabilizing results in another recurrent 
sandpile; so all sandpiles in C v are recurrent. Note that 

{V + Ps)° = (rj + e s + /3 S )° = r) ® j s 

where j s = ( e s + f3 s )° is the recurrent representative of (3 S . It follows that 

Cl, = r] © (7s) 

where ( 7 ^) denotes the cyclic subgroup of Rec(G, s) generated by 7 *. 

To investigate these cosets of ( 7 ^), we recall the period vector introduced by Bjorner 
and Lovasz. 

Definition 2.6. [ ] Given a graph G with total Laplacian A, a vector p £ N n is 
called a period vector for G if p 7 ^ 0 and Ap = 0. A period vector is primitive if 
the greatest common divisor of its entries is 1 . 

In other words, a period vector p has the property that firing each vertex v £ V a 
total of p(u) times results in no net movement of chips. The following lemma sums 
up some useful properties of period vectors. 

Lemma 2.7. [ , Prop. 4.1] A strongly connected multigraph G has a unique primitive 
period vector 7 tq- All entries of hq are strictly positive, and all period vectors of G 
are of the form kite for k = 1,2,.... Moreover, if G is Eulerian, then ttq = 1 . 

A consequence of the strict positivity of 7 tq that we will use several times is that 
AZ n = AN n . 

We now introduce a very special period vector. Recall that n(v) denotes the number 
of spanning trees of G oriented toward v. 
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Lemma 2.8 ([1, 9]). A k = 0. 

Recall the Pham index M = Mq, defined as the greatest common divisor of the 
spanning tree counts {k(u)|'u G R}. By Lemmas 2.7 and 2.8, the vector n = jjk is 
the unique primitive period vector of G. 

Next we argue that n(s) = ord( 7 s ), the order of g s in the group Rec(G, s). Fixing 
a positive integer m, we have that m/3 s is trivial in K(G, s ) if and only if there exists 
y G Z n_1 such that m/3 s = A s y. Setting x = (y, 0) G Z n and noting that (3 S is the 
restriction of —A 5 S to the nonsink vertices, such y exists if and only if there is a 
vector x G Z n such that x(s) = 0 and A(x + mS s ) = 0. (The equality in the sink 
coordinate of the first equation follows from the equality in the nonsink coordinates 
because the sum of all the coordinates is 0.) Setting p = x + mS s , such x exists if 
and only if there is a vector p G Z n such that 

Ap = 0 and p(s) = m. 

By Lemma 2.8, since the kernel of A is one-dimensional and m is positive, such p 
must be a positive integer multiple of the primitive period vector n. Such p exists 
if and only if 7 r(s) divides m. Thus 7 r(s) is the order of /3 S in I\(G,s), which by 
Theorem 2.5 is the order of in Rec(G, s). Recalling that 7 r(s) = k(s)/M, we 
conclude the following. 

Lemma 2.9. [22, Lemma 6 ] For any choice of sink s, we have that 

ord( 7 s ) = k(s)/M = #Rec (G,s)/M 

Thus, M = ffRec(G,s)/ ( 7 *) is the number of distinct cosets of ( 7 ^) in Rec(G, s). 

2.3. Comparison of sandpile groups with and without sink. We now investi¬ 
gate the structure of the quotient group Rec(G, s)/ ('y s )- Recall that (3 S is the sandpile 
(3 s (v) = d sv , where s is the designated sink vertex and that /3 S is the equivalence class 
of j3 s in K(G, s ). As before we write Zg for the group of vectors in Z n with coordinates 
summing to 0 . 

Theorem 2.10. For any strongly connected multigraph G and any vertex s, 

Rec(G, s)/ < 7s ) =* K(G, s)/ {J s ) ^ Z^/AZT 

The meat of the proof for this theorem is packaged in the following workhorse 
lemma. To translate between sandpiles and total configurations, we introduce some 
notation: Given a vector x G Z", we denote by x the restriction of x to the nonsink 
vertices; and given g G Z n_1 , we write f° r the extension of g to Z n such that 
\Vk\ = k. 

Lemma 2.11. Let ct,t gZ" with |cr| = |r|. Then the following are equivalent. 

(1) a = r mod AZ n 

(2) <7 = f mod A s Z n_1 + Z/3 S 

(3) (a + e s )° G (f + e s )° ® (g s ) 
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Proof. (1 2) Assume (1), and let m = \a\ = |r|. Recall that o^ denotes the 

extension of o to Z n such that |cjfc| = k. We observe that (1) holds if and only if 
there is an x E Z" 1 such that o = r — Ax. If cr = r — Ax, then 


a = t — Ax = t — A 


x 

0 


-A 


0 

x(s) 


= T — 


A s x 

a 


- x(s) c s 


where c s denotes the column of A corresponding with the sink and a is the dot product 
of the nth row of A with (x, 0). Since /3 s (z) = —c s (i) for each i s, it follows that 
o = t — A s x + x(s)/3 s . Going the other way, we assume that d = t — k(3 s — A s x for 
some ken and x E N n_1 . Let o' be the total configuration 


Then o'(%) = o{i) for all i ^ s and \o'\ = |r|. Since cr(s) is determined by |oj and 
|cr| = |r|, we have that o' = o. 

(2 -<=> 3) Note that (3) is equivalent to the existence of an x such that 
(■o + e s )° = (f + e s )° + x(s) (/3 S + e s )° - A s x 
which in turn is equivalent to the congruence o = t mod A s Z n_1 + Z f3 s . □ 


Proof of Theorem 2.10. Define a map (f : K(G,s)/ (/ 3 S ) —> sending 

fj mod (/3 S ) i-a r ]o mod AZ n . 

Let E Z n_1 . If fj = ^ mod (/3 S ), then by Lemma 2.11 we have that ?/o = 
mod AZ n , so that cf is well-defined. The equation % + Co = iv + £)o immediate 
from the definition, so that (f is a homomorphism. The map is also surjective, 
since for each o E Zq there is a corresponding o E Z n_1 , and ( o mod (/3 S )) = 
o mod AZ n . We now show that cj) is injective to complete the proof that cj) is an 
isomorphism. Suppose that o = r mod AZ n . Then by Lemma 2.11 we have that 
o = t mod and the theorem is proved. □ 


2.4. Eulerian and CoEulerian Graphs. We now prove the two results stated in 
the introduction. We also supplement Proposition 1.1 with three equivalent condi¬ 
tions about the sandpile group. 


Proposition 2.12. The following are equivalent for a strongly connected directed 
multigraph G = (V,E). 

(1) ker(A : Z y -»■ Z v ) = Zl. 

(2) Mq = k{s) for all sEh 

(3) G is Eulerian. 

(4) For all s E V, the element /3 S is trivial in the sandpile group K(G,s). 

(5) K(G,s ) “ Zq "/ATX for all s E V. 

(6) K(G, s) = K(G , s') for all s, s' E V. 

Proof. (1 3) We have A1 = 0 if and only if the indegree of each vertex equals 

its outdegree. By [25, Theorem 5.6.1], for G strongly connected this degree condition 
is equivalent to the existence of an Eulerian tour. 
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(1 =>• 4) By definition, j3 s is the restriction of —A 5 S to the nonsink coordinates. 
If A 1 = 0 then -A S s = A(jj v ^ s S v ), so f3 s = A S (X^ S € A s Z n_1 . 

(4 => 5) This follows directly from Theorem 2.10. 

(5 ==> 6 ) Trivial. 

(6 =>• 2) If K(G,s ) = K(G,s'), then equating orders yields k(s) = k(s'). If this 
holds for all vertices s and s', then Mq := gcd{K(u)|u € V} = k(s) for all s £ V. 

(2 => 1) If all coordinates of k are equal, then Al = 0 by Lemma 2.8. □ 

In particular, the sandpile group K(G,s ) is independent of the choice of sink if 
and only if G is Eulerian (the “if” direction is well known [ 7, Lemma 4.12]). 

Proof of Theorem 1.2. (3 =>■ 4) We prove the contrapositive. Assume there is 

a sink s such that K(G,s ) ^ (/3 S ), and fix the number of chips on G to be m = 
#E — ffV. Our assumption implies that there are two distinct cosets C\ and C 2 of 
( 7 S ) such that <r max £ C\. Choosing an 77 £ C 2 , we remark that r] m stabilizes if and 
only if rj m = cr max mod AZ n (since cr max is the only stable total configuration with 
m chips). By Lemma 2.11, this congruence holds if and only if i] £ cr max © ( 7 ^) = Ci, 
so we see that 7j m does not stabilize. 

(4 => 5) Trivial. 

(5 => 3) Let a be a total configuration with |cr| < ffE — ffV. We first write 
a = t — 6 for some r, 5 £ Z n where |r| = ffE — ffV and 6 > 0. Now if (f3 s ) = K(G, s ) 
for some vertex s, then f = a max mod AsZ”^ 1 + Z j3 s so that r = cr max mod AZ n by 
Lemma 2.11. It follows that a = cr max — 5 mod AZ n . Using that AZ n = AN n , we 
conclude from Lemma 2.1 that a stabilizes. 

(1 5) This equivalence follows from Theorem 2.10. 

(2 5) This equivalence follows from Lemma 2.9. □ 

2.5. Graphs that are both Eulerian and coEulerian. We conclude this section 
by characterizing the graphs that are both Eulerian and coEulerian. A strongly 
connected graph G without loops is called a directed cactus [28, 23] if each edge 
of G is contained in a unique simple directed cycle. Let us call this the “unique 
cycle property” (UCP). As the proof of the next Proposition will show, the UCP is 
equivalent to the following “unique path property” (UPP): For any pair of vertices 

x, y £ V there is a unique simple directed path in G from x to y. (A simple cycle or 
path is one with no repeated vertices; in particular, a 2 -cycle consisting of an edge 
and its reversal is simple.) 

Directed cacti are in some sense analogous to trees. In particular, a bidirected 
graph has the UCP if and only if it is a bidirected tree. 

Proposition 2.13. Let G be a strongly connected finite graph without loops. Then 
G is both Eulerian and coEulerian if and only if G is a directed cactus. 

Proof. Supposing that G is both Eulerian and coEulerian, we have n(y) = 1 for 
all y £ V. By Wilson’s algorithm [27], any oriented spanning forest, one of whose 
components is oriented toward y, can be completed to a spanning tree oriented toward 

y. Given x, y £ V and a simple path P from x to y, completing P in this manner 
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results in the unique spanning tree T y oriented toward y. Therefore all simple paths 
from x to y are contained in T y , so G has the UPP. 

Next observe that for each directed edge e = (y, x ) there is a bijection between 
simple directed paths P from x to y and simple directed cycles P U {e} containing e. 
Hence the UPP implies the UCP. 

Finally, supposing that G has the UCP, we will compete the proof by showing that 
for each vertex y there is a unique spanning tree of G oriented toward y, so that G 
is both Eulerian and coEulerian. Let C be the set of simple directed cycles in G, and 
consider the undirected bipartite graph T on vertex set V U 6 whose edges are the 
pairs {v,C} such that vertex v lies on cycle C. The UCP implies that T is a tree. 
Now we can manifestly describe the unique spanning tree of G oriented toward y. 
Namely, for each cycle C let ec be the outgoing edge from x in C , where (C, x,.. ., y) 
is the unique path from C to y in T. The remaining edges E — {ec\C € C} form 
a spanning tree of G oriented toward y. Moreover any such spanning tree T must 
contain all edges of G of the form (x, x') where (C, x, C', ..., y) is the path from any 
cycle C to y in 7 and x' € C 1 , else there would be no path from x to y in T. By the 
UCP, edges ec and (x, x') are distinct since they belong to distinct cycles. Therefore 
T must omit all of the edges ec, and hence T is unique. □ 


3. Computational complexity 

Bjorner and Lovasz [5, Corollary 4.9] showed that the halting problem for chip-firing 
can be decided in polynomial time for simple Eulerian graphs. By Theorem 1.2, it 
can be decided in linear time for coEulerian multigraphs. The purpose of this section 
is to show that despite these two easy cases, the problem is NP-complete for general 
directed multigraphs. 

To see that it is in NP, let a be a nonnegative halting chip configuration on a 
strongly connected directed multigraph G = (V,E), and let x(u) be the number of 
times vertex v fires. By Lemma 2.1 the vector x is a certificate that a halts. Why 
does this certificate have polynomial size? By Lemma 1.4 we have x(u) < n(v) for 
some vertex v. Moreover for any directed edge (ui,U 2 ) the vertex U 2 receives at least 
x(ui) chips from u\ and so U 2 fires at least x.(ui)/d U2 times. For any vertex u , by 
inducting along a path from u to v we find that x(«) < Dx(r) where D = n^ev d w 
is the product of all outdegrees. By Lemmas 2.7 and 2.8 relating the primitive period 
vector 7r to the spanning tree count vector k. we have n(v) < k(v) < D, so all 
entries of x are at most D 2 . Noting that log D < v£V log d uv , which is the size of 
description of the adjacency matrix, we conclude that the halting problem for 
chip-firing is in NP. 

To show that it is also NP-hard, our starting point is the following decision problem 
considered by Amini and Manjunath [ ]. 

Nonnegative rank: 

Given a basis of an (n — l)-dimensional lattice L C Zq and a vector a € Z n , 

Decide whether there is a vector t € N n such that a — r £ L. 
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If there exists such a r, then a is said to have nonnegative rank relative to L. In 
[2, Theorem 7.2] nonnegative rank is shown to be NP-hard by reducing from the 
problem of deciding whether a given simplex with rational vertices contains an integer 
point. (To give a little context, the term “rank” is inspired by the Riemann-Roch 
theorem of Baker and Norine [ ]. Asadi and Backman [3] extend parts of the Baker- 
Norine theory to directed graphs. Kiss and Tothmeresz [ ] show that computing 

the Baker-Norine rank—a harder problem than deciding whether it is nonnegative—is 
already NP-hard when L is the Laplacian lattice of a simple undirected graph.) 

The link between chip-firing and nonnegative rank is provided by the following 
variant of a theorem of Perkinson, Perlman and Wilmes [20]. 

Theorem 3.1. Given an (n — l)-dimensional lattice L C Zq, there exists a strongly 
connected multigraph with Laplacian A such that 

L = AZ n . 

Moreover, A can be computed from a basis of L in polynomial time, and all entries 
of A are bounded in absolute value by nd where d = detL. 

The inspiration for Theorem 3.1 is [20, Theorem 4.11], which expresses an arbitrary 
(n — l)-dimensional lattice in Z n_1 as a reduced Laplacian lattice A s Z n_1 . Modifying 
its proof to express L C Zq as a total Laplacian lattice is straightforward; we give 
the details below. 

In our application it will be essential to compute the Laplacian matrix A from a 
basis of L in polynomial time (in the length of description of the basis). It is not 
evident whether [ , Algorithm 4.13] runs in polynomial time, due to possible blow up 

of the matrix entries in repeated applications of the Euclidean algorithm [ 1,1 ]. As 
detailed below, this numerical blow up can be avoided by the usual trick of computing 
modulo the determinant d. 

To see how we will apply Theorem 3.1, note that strong connectivity implies 

AZ n = AN n 

since the period vector of Lemma 2.7 is strictly positive. Thus, a vector a G Z n has 
nonnegative rank relative to L = AZ n if and only if there exists x € N n such that 

a + Ax > 0. 

Now by Lemma 2.1, such an x exists if and only if the chip configuration cr max — a 
stabilizes. To summarize, a polynomial time computation of A given a basis for L 
yields a polynomial time Karp reduction from NONNEGATIVE RANK to the HALTING 
problem FOR chip-firing on a finite directed multigraph, showing that the latter 
is NP-hard. 

Corollary 3.2. the halting problem FOR chip-firing is HP-complete. 

It remains to prove Theorem 3.1. Recall that an m x m integer matrix U is called 
unimodular if det U = ±1. Any nonsingular square integer matrix A has a Hermite 
normal form 


H = AU 
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where U is a unimodular integer matrix, and H = ( hij ) is lower-triangular with 
integer entries satisfying 


0 < ha, 1 < i < m 

0 < < ha, 1 < j < i < m. 

The existence and uniqueness of H was proved by Hermite [16]. The Hermite normal 
form is useful to us because H can be computed from A in polynomial time [12] and 
HZ m = A(U7/ n ) = AZ m by the unimodularity of U. Let 

m 

d = | det A\ = det H = JJ ha. 

2=1 

We will use the following observations about the column span AZ m . 

Lemma 3.3. [ 12, Cor. 2.3] Let B be a lower triangular m x m matrix whose columns 
are in and whose diagonal entries satisfy bn = ha for all i. Then L>Z m = HZ m . 


Lemma 3.4. [12, Prop. 2.5] dZ m C AT, m . 

We will apply these lemmas with rn = n — 1. Note that an n x n integer matrix 
is the total Laplacian of a directed multigraph if and only if (i) the entries of each 
column sum to zero, (ii) the diagonal entries are nonnegative, and (iii) the off-diagonal 
entries are nonpositive. 

Given an n x (n — 1) integer matrix whose columns are a Z-basis of the (n — 1)- 
dimensional lattice L C Zq, let A be the result of removing the last row of M. Since 
each column of M sums to zero, A is nonsingular. Let H = AU be the Hermite 
normal form of the (to — 1) x (to — 1) nonsingular matrix A. 

The hypotheses of Lemma 3.3 are trivially satisfied when B = H; and if B satisfies 
the hypotheses of Lemma 3.3 then by Lemma 3.4 it will continue to do so if we 
subtract d from an entry below the diagonal. Using this operation we can make the 
entries immediately below the diagonal sufficiently negative so that the sum of the 
entries in each column is nonpositive. Namely, let B = (b t j) where 


bij 


h^ — kjd, i = ,7 + 1 
hij , else 


and kj for each j = 1,..., to — 2 is a nonnegative integer such that 

m 

(kj — 1 )d < ^2 hij < kjd. 

2=1 

Now let 


-\-d 

—/in 

0 

0 

0 

0 ■ 

0 

+ 

— /&22 

0 

0 

0 

0 

— 

+ 

-/l 33 

0 

0 

0 

— 

— 

+ 

—L44 

0 

0 

_ 

_ 

_ 

_ 

fomm 

. -d 

- 

— 

— 

— 

-\~hmm _ 


( 2 ) 
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be the nxn matrix with upper right corner —B, the column vector dei~de n appended 
on the left, and a row appended on the bottom such that the entries of each column 
sum to zero. By the choice of kj in (2), the bottom row of A is nonpositive, except 
for its rightmost entry h rnm . Therefore A satisfies the conditions (i)-(iii) above. Since 
the entries immediately above the diagonal of A are negative, as is A„i, the matrix 
A is the Laplacian of a strongly connected multigraph (it has the Hamiltonian cycle 
1 —> n —> n — 1—»■■■■—>-1). By Lemma 3.4 the first column of A belongs to L. 

Moreover, since both L and AZ n are contained in Zq and BZ n_1 = HZ’ 1-1 , the 
integer span of the remaining columns of A is L. Thus L = AZ n . Since each entry 
of H is at most d = Yl ha, each entry of A has magnitude at most nd, completing 
the proof of Theorem 3.1. 

3.1. Simple directed graphs. Let us point out a sense in which the NP-hardness of 
Corollary 3.2 is rather weak: the directed graphs for which the HALTING PROBLEM FOR CHIP-FIRING 
is hard may have large edge multiplicities. This is because the Laplacian A of The¬ 
orem 3.1 may have large entries, which in turn is because the lattice L in a hard 
instance of nonnegative rank has large determinant. An interesting question is 
whether the halting problem for chip-firing remains NP-hard when restricted 
to simple directed graphs, those with edge multiplicities in {0,1}. Does the hardness 
arise from directedness or from large edge multiplicities (or both)? The following 
table summarizes what is known. 



simple graphs 

multigraphs 

coEulerian 

P (Theorem 1.2) 

P (Theorem 1.2) 

bidirected 

P (Tardos 26 ) 

? 

Eulerian 

P (Bjorner-Lovasz [5]) 

? 

strongly connected 

? 

NP-complete (Cor. 3.2) 


Table 1. Complexity of the halting problem for chip-firing for 
eight different classes of strongly connected directed graphs. 
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